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Abstract. This note concerns the relationship between conditions on cost functions and 
domains and the convexity properties of potentials in optimal transportation and the conti- 
nuity of the associated optimal mappings. In particular, we prove that if the cost function 
satisfies the condition (A3), introduced in our previous work with Xinan Ma, the densities 
and their reciprocals are bounded and the target domain is convex with respect to the 
cost function, then the potential is continuously differentiable and its dual potential strictly 
concave with respect to the cost function. Our result extends, by different and more direct 
proof, similar results of Loeper proved by approximation from our earlier work on global 
regularity. 



We continue our investigation on the regularity of potential functions in the optimal 
transportation problem [MTW, TW]. In this paper we prove strict convexity and 
regularity of potential functions for non-smooth densities. The strict convexity and 
regularity for solutions to the Monge-Ampere equation were established by Caffarelli 
[CI]. For the reflector design problem, which is a special optimal transportation problem 
[W2], these results were obtained in [CGH]. For more general optimal transportation 
problems, the C^'" regularity has been obtained by Loeper [L]. In this paper we prove 
the strict convexity for potential functions, and obtain the regularity under weaker 
conditions on the domains. We also use our results to plug a gap in [MTW] pertaining 
to the use of a comparison argument in the proof of interior regularity. 

Let 0,0* be two bounded domains in M", and f,g be two nonnegative integrable 
functions on O, O* satisfying the mass balance condition 



Let (u, v) be potential functions to the optimal transportation problem, namely (w, v) is 
a maximizer of 



1. Introduction 



(1.1) 





sup{/((^,V') : e K}, 
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where 



(1.3) IM)= j f{xMx)+ f ^(yMy), 

Jci Jq* 

K = {{^, ^) e C0(1]) X C\Q*) : ^(x) + i^iy) < cix, y)}. 

We assume that the cost function is smooth, c G C°°(M'^ x R"^), and satisfies (Al)-(A3) 
below. It is known [C4, GM] that there is a maximizer {u,v) to (1.2), which is also 
unique up to a constant if /, g are positive. The potentials (w, v) are semi-concave and 
satisfy 

(1.4) u{x)= M {c{x,y)-v{y)}, 

yen* 

v{y) = inf {c(a;,y) -^(a;)}. 

The optimal mapping T : x & Q ^ y & can be determined, a.e. in Q, by 

(1.5) Du{x) = Da;c{x,y), 

where y also attains the infimum in (1.4). If tt e C^(ri), it satisfies the equation 

f 

(1.6,) det(i:>^c- i:'^^) = |detcij|-^ in 

where Cj = d^iC, Cij — ^x^^y.c. The main result of this note is 

Theorem 1. Suppose that O* is c- convex with respect to Q, that the cost function c 
satisfies (Al)-(A3), and that f,g satisfy 

(1.7) Ci</,^<C2 

for some positive constants Ci, C2. Then v is strictly c-concave and u is smooth. 

We refer the reader to §2.1 for definitions of various convexity notions relative to cost 
functions. By approximation and the uniqueness of potential functions (when f,g > 0), 
condition (1.7) can be weakened to < < C and f,g > 0. From the smoothness 
of u it follows that the optimal mapping T is continuous. 

For the reflector design problem (in the far field case). Theorem 1 was obtained in 
[CGH]. For the optimal transportation problem, the smoothness of u essentially be- 
longs to Loeper [L]. His proof uses approximation by globally smooth solutions to the 
optimal transportation problem, established in [TW] , and accordingly assumes additional 
conditions such as the c-convexity of both domains fl and fl*. See also remarks after 
Corollary 1 in §2.8. Here we give a more direct proof. Our proof is also completely differ- 
ent from that in [CI], in which Caffarelli proved the strict convexity and regularity for 
solutions to the Monge-Ampere equation with constant boundary condition. For higher 
regularity, the interior and global C^'" estimates for solutions of (1.6), in the case of the 
quadratic cost function, were established in [C2, C3, Ul], and earlier in [D] for n = 2. 
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For cost functions satisfying (Al)-(A3) below, the regularity of potential functions was 
obtained in [MTW, TW]. 

By approximation and the regularity in Theorem 1, it follows that the potential 
function u is fully c-concave if /, (7 > and c satisfies (Al)-(A3). In particular, every 
local c-support is a global one. See §2.8 for more discussion. 

The assumption (Al)-(A3) follows. 

(Al) For any z E M", there exists a unique y = y{x, z) such that Dxc{x, y) — z. 

(A2) For any x, y G det{cij(a:, y)} ^ 0. 

(A3) For any x eVL,y eVt* , and ^, ry e M'^ with ^ -L r/, 

^-^i,3,k,l,p,q,r,s 

where cq is a positive constant, and (c*'-') is the inverse matrix of {cij). 

We also assume the above conditions hold after exchanging x and y. Formula (1.8) is 
equivalent to 

(1-9) d^kZ,Cij{x,y)^iCj'nkVi < -colCI^I^?!^ 

for all a; G G O* and ^ _L G M"^, where y = y{x^z) is given by (Al), which is 
smooth in x and z by (A2). Note that for global regularity in [TW] and the subsequent 
application in [L], condition (A3) can be relaxed to its degenerate form, cq = 0, called 
A3w in [TW]. 

We divide the proof of Theorem 1 into several short sections. We first introduce in §2.1 
various notions of convexities and concavities relative to the cost function c. We then 
indicate in §2.2 a geometric property of (A3) (see also [L]). In §2.3 we give an analytic 
formulation of the c-convexity of domains. A geometric characterization of c-convex 
domains (under condition A3) is given in §2.4. In §2.5 we prove that a local c-concave 
function is fully c-concave if the domain is c-convex and the cost function c satisfies (A3). 
This is a key ingredient in the proof of Theorem 1. In §2.6 we show that u is if and 
only if V is strictly c-concave. We then prove Theorem 1 in §2.7. Various remarks are 
given in §2.8. 



2. Proof of Theorem 1 

2.1. Convexities relative to cost functions [MTW]. Let u be a semi-concave func- 
tion in f2, namely u — C\x\'^ is concave for a large positive constant C. The supergradient 
d^u [GM] and c-supergradient d'^u are defined by 

(2.1) 

d'^u{xQ) — {p : u{x) < u{xo) + p ■ {x — xq) + o{\x — xq])}, 

(2.2) 

d^u{xo) — {y eM.^ : u{x) < c{x, y) — c(a;o, y) + u{xq) -\- o{\x — xo\)} 
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for X near xq, where xq G ft. For a set e f2, we denote d^u{E) = \Jx£Ed'^u{x) and 
d+u{E) = Ua;eEd+u{E). By (1.5) we have 

(2.3) Cx{xo, d^u{xo)) = d^u{xo). 

Note that d'^u{xo) is a dosed, convex set. Hence d^u{xo) is cfosed and c-convex with 
respect to a^o- 

We may extend the above mappings to boundary points. Let xo € dfl be a boundary 
point, wc denote d~^u{xQ) = {p G M"" : p = hmfc^oop^}, where pk G d~^u{xk) and {.x^} 
is a sequence of interior points of Q such that xq, and let 9+u(a;o) be given by 

(2.3) . 

The c-normal mapping is defined by 

(2.4) Tu{xo) = {y e R" : u{x) < c{x, y) - c{xo, y) + u{xo) for all x e Q}. 

Note that Tu{xo) C d^u{xo). 

c-support: Let w be a semi-concave function in Q. A local c-support of u at G ^ is 
a function of the form 

h = c(-,yo) + oo, 

where ciq is a constant and yo G M", such that u{xo) = h{xQ) and u(x) < h{x) near xq. 
If u{x) < h{x) for all x & Q, then /i is a global c-support (or c-support for short) of u at 
xq. If /j. is a local c-support of u at a;o, then yo G 5+tt(a;o) and ao = ^(aJo) — c(a;o) Ho)- 

c-concavity of functions: We say a semi-concave function u is locally c-concave if for 
any point xq & 0, and any j/ G d^u{xo), h = c(-, y) — c{xo, y) -\- u{xq) is a local c-support 
of tt at xq. We say u is c-concave if for any point xq G fi, there exists a global c-support 
at Xq in O. We say u is strongly c-concave if it is both locally c-concave and c-concave. 
We say u is fully c-concave if it is locally c-concave and every local c-support of u is a 
global c-support. We say u is strictly c-concave if it is fully c-concave and every c-support 
of u contacts its graph at one point only. 

c-segment: A set of points £ C R"' is a c-segment with respect to a point yo G if 
Dyc{£, yo) is a line segment in R". 

c-convexity of domains: We say a set U is c-convex with respect to another set V if 
the image Cy{U,y) is convex for each y eV. Equivalently, U is c-convex with respect to 
V if for any two points xq, Xi E U and any y E V, the c-segment relative to y connecting 
Xo and Xi lies in O. By (2.3), a c-convex domain is topologically a ball. 

By definition, a c-concave function u can be represented as [GM] 

u{x) = inf{c(x, y) - a{y) : y G T„(f2)}, 

where a is a function of y only. By (1.4), a potential function u is c-concave [GM]. Our 
Theorem 1 implies that u is furthermore fully c-concave under assumption (A3). If u is 
C^, then local c-support is unique and c-concavity is equivalent to full c-concavity. We 
also remark that a potential function may not be locally c-concave in general. 
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Similarly we can define c*-segment, c*-support, c*-convexity and c*-concavity by ex- 
changing variables x and y [MTW] . In this paper, we will generally omit the superscript 
* when the meaning is clear. 

2.2. A geometric property of (A3). Let 2/0,2/1 be two points in O*. Let yoyi be 
the c-segment relative to a point xq G Q, connecting 2/0 and 2/1- By definition, 

(2.5) y^i = {yt : c^{xo,yt) = Pu te [0, 1]}, 
where pt = tpi + (1 - t)po, and po = Cx{xo, 2/0), Pi = Cx{xo, 2/1). Let 

(2.6) hi{x) = c{x,yi) - tti, i = 0,l, 

ht{x) = c{x,yt) - at, 

where t G (0, 1), ao, ai and at are constants such that ho{xo) = hi{xo) = ht{xo). Suppose 
(A3) holds. Then for x xq, near xq, we have the inequality 

(2.7) ht{x) > mm{ho{x),hi{x)}, 

which is crucial for the remaining analysis of this paper. 

Inequality (2.7) follows from (1.9). Indeed, by a rotation of axes, we assume that 
Pi = Po + ^^n, where 5 is a positive constant and Cn = (0, • • • , 0, 1) is the unit vector in 
the x^-axis. By (1.9), 

(2.8) -^Cij{x,y{x,pt))iiij < -co5^ 

for any unit vector ^ orthogonal to e^, where y = y{xo,pt) is given in (Al). Now (2.7) 
follows from (2.8); for details see [L]. 

2.3. An analytic formulation of the c-convexity of domains. If O is c-convex 
with respect to fi*, by definition, Cy{0., y) is convex for any y . Suppose G dVt and 
locally 30. is given by 

(2.9) Xn = p{x') 

with Dp{0) = such that e„ is the inner normal at 0, where x' — {xi, ■ ■ ■ ,x„_i). Then 
at (0, y) (for a fixed point y), 

9xiCy = Ci,y + Cn,yPxii 

(2.10) dx^xj {Cy,l) = {Cij,y, 7) + {Cn,y,l)PxiXj > 0, 

where 7 is the inner normal of Cy{Q,y) at Cy{0,y). We may write (2.10) explicitly, 
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Make the linear transformation 

Vk = akiyi- 

Then we have 

(2-11) Cij^y^akia;;^]^rn + Cn,U(^kiamnniPij > 0. 

Let Qij = Cij{0,y). Then Cxi,yj — Sij and 7 = Cn- We obtain 
(2-12) Cy,^„ +pxixj > 0, 

which is equivalent to 

(2.13) c,,-,,c'='- + p,,,. >0. 

Let if G C^(0) be a defining function of O. That is (/? = 0, |V<^| on dO, and (p < 
in Q. From (2.13), we obtain an analytic formulation of the c-convexity of Q relative to 

n* [TW], 

(2.14) - c'''\x,y)cij,k{x,y)(piix)] > ^ x edn,y en*. 

Conversely, if O is simply connected and (2.14) holds, then O is c-convex. Following 
[TW], we call fl uniformly c-convex with respect to Q* if the matrix in (2.14) is uniformly 
positive. 

2.4. Geometric properties of the c-convexity of domains. Let j/cJ/i be any 
two given points in Q*. Denote 

(2.15) Af = Myo,y,,a = {xeW : c{x, yo) = c{x, yi) + a}, 

where a is a constant. Assume that the origin e M and locally is represented as 

(2.16) Xn = v{x') 

such that 77(0) = and Dr]{0) = (obviously rj also depends on yo,yi and a). Then 

(2.17) c{x',r]{x'),yo) = c{x' ,r]{x'),yi) + a. 
Differentiating (2.17) gives 

Cj(0, yo) + CnVi = Cj(0, yi) + CnVi, 

(2.18) Cjj(0,yo) + CnVij = Cy(0,yi) + CnVij- 

We obtain 

(2.19) [cn(0,yi) -Cn(0,yo)]%- + [cij{0,yi) -Cij(0,yo)] = 0. 
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Now let yo be fixed but let yi and a vary in such a way that yi — > yo and the set 
Af = Myo,y-^,a, given by (2.15), is tangential to {xn = 0}, namely ?7(0) = and Dr]{0) = 0. 
By a linear transform as in §2.3 we assume that Cxi^y^ = Sij at a; = and y = yo- Then 



yi - yo 

\yi - yo\ 



We obtain 



(2.20) Cn,n(0,yo)r/°- + Cy, n(0,yo) = 0, 

where rj^ is the limit of rjy^^y^^a as yi and a vary as above. 

Now let O be c-convex, (uniformly c-convex), with respect to O*. Suppose that 90 
is given by (2.9) with -Dp(O) = so that dVL is tangential to N at the origin. Then by 
(2.12)and (2.20) we obtain 

(2.21) D'^{p-rf)>{), (>0), atO. 



From (2.21) we obtain some useful geometric properties of c-convex domains, assuming 
that c satisfies (A3). 

First, if Vt is c-convex with respect to O*, then for any compact subset G C O*, O 
is uniformly c-convex with respect to G. Indeed, let yo^yi be two points in O*. Let 
Pq = Cx{0,yo), pi = Cx{0,yi). For t e (0, 1), let pt = tpi + (1 - t)po and yt satisfy 
Ca;(0, yt) = Pt- Then the set {pt : < t < 1} is a line segment and the set {yt : < t < 1} 
is a c-segment. Suppose pi = Po + Sen for some d > 0. Let ht{x) = c{x,yt) + at, where 
at = c(0,yo) - c(0,yt) such that ht{0) = ho{0) for all t e [0,1]. Let Aft = {x e R"" : 
ht{x) — ho{x)} and suppose that near 0, J\ft is given by Xn = r]t{x'). Since Ca;(0, yt) — pt, 
we see that Aft is tangential to {xn = 0}, namely D{rit — ?7t')(0) = 0- By (2.8) we have 
furthermore the monotonicity formula 

(2.22) D\r]t- m'){0)>0 



for any t > t' and t, t' G [0, 1]. Geometrically it implies that Aft lies above Aff ii t > t', 
namely rjt{x') > rjt'{x') for x' near 0, and equality holds only at x' = 0. Consequently 
we obtain from (2.21) the strict inequality 

(2.23) D^{p-r])>0 at 0. 



From (2.23) we obtain the above mentioned property. 

Next, for any j/O; J/i £ if (given in (2.15)) is tangent to dfl at some point xq 
and if Q is c-convex with respect to O*, then the whole domain Q lies on one side of Af. 
Indeed, we may assume xq = 0, locally dfl and Af are given respectively by (2.9) and 
(2.16), such that Dp{0) = Dr]{0) = 0. By (2.23), p{x) > r]{x) for x near 0, x ^ 0. Denote 



U = Uyo,yi,a = {a; e R" : c{x, yo) < c{x, yi) + a} 
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so that M = dU. If fl does not lie on one side of jV, namely Cl is not contained in C/, then 
Q — U contains two disconnected components (one is the origin). Since 90 is a closed, 
compact hypersurface, we decrease the constant a (shrinking the set U) until a moment 
when two components of fl — U meet each other at some point x* G dfl. But since J\f is 
tangent to dfl at x*, we reach a contradiction by (2.23). 

It follows that if O is c-convex with respect to O*, then 
(2.24) ^ = C]Uyo,yua, 

where the intersection is for all yo,yi £ O* and constant a such that Uy^^y^^a ^ O. 
However, we don't know if the converse is true, namely whether Q. is c-convex with 
respect to Q,* if it is given by (2.24). 

The above properties also extend to cost functions satisfying A3w and uniformly c- 
convex domains. 



2.5. Local c-support is global. Let u be a locally c-concave function in Vt with 
9+i(.(0) C O*. Suppose Q, is c-convex with respect to O*. Let h{x) = ha{x) — c{x, yo) + a 
be a local c-support of u ai, xq. Then /i is a global c-support of u, namely 

(2.25) u{x) < h{x) V X e O. 

Indeed, if this is not true, then for £ > small, the set {a; G O : u{x) > ha-ei^)} 
contains at least two disconnected components. We increase e (moving the graph of h 
vertically downwards) until at a moment e = Eq > 0, two components first time touch 
each other at some point xq. If xq is an interior point of Q, by definition ha-so cannot be 
a local c-support at xq, which implies that yo ^ d^u{xo). We claim that for a sufficiently 
small r > 0, j/o does not lie in the set d^u{Br{xo)) either. Indeed, if hk = c(-, j/o) + cbk 
for k = 1,2, ■ ■ ■ is a sequence of local c-support of u at x^ and if x^, — > xq, oq, we have 
yo ^ dt^i^o) ^ semi-concave. Hence = c(-,yo) + ao is a local c-support of u at 
Xo- This is a contradiction. 

Therefore ha-eo ^ transversal near xq, and for e < Eq, close to Eq, locally the 
set {x e O : u{x) > ha-s{x)} cannot contain two disconnected components. Hence xq 
must be a boundary point of Q. 

In case xq is a boundary point of Q, we will also reach a contradiction by (2.23). 
Similarly as above, ha-eo cannot be a local c-support at xq, namely yo ^ d^u{xo). 
Without loss of generality let us assume that xq — and locally dfl is tangent to {xn = 0} 
such that e„ is an inner normal of O at 0. As before we also assume that Cj^^^y^ — Sij at 
X = and y = yo. Let po = Cx{0,yo)- By subtracting a linear function of x from both 
c and u, we assume that {xn = 0} is a tangent plane of h at 0. Then we have po — 0, 
u{0) = 0, and u{x) < o(|a;|) as x — > 0. Since yo ^ d^u{xo), we have 

(2.26) P = lim ju{ten) < 0. 

Let pi ~ (3en, Pt = tpi + (1 ~ i)Po ^ot t G [0, 1], and let yt G be determined by 
Ca;(0,2/t) = Pt- Then pi G d^u{Qi) and yi G d'^u{Q). Since u is locally c-concave. 



c{x, yi) + 6 is a local c-support of u at 0, where 6 is a constant such that c(0, yi) + b = 
u{0) = 0. Denote 

U = {xeW: ha-soix) > c{x, yi) + b} 

and J\f = dU = {x eW^ : ha-eo {x) = c{x, yi) + h}, such that e jV. Since Q, is c-convex, 
by (2.23) we see that Q, <zU. But since c(a;, yi) + 6 is a local c-support of u at 0, we have 

(2.27) h 

for X e n, near the origin. We reach a contradiction by our choice of Sq. This completes 
the proof of (2.25). 

From the above proof, we see that if /t is a c-support of then the contact set 
{a; e f2 : h{x) = u{x)} cannot contain two disconnected components (or points). In 
other words, the contact set is connected. 

We also remark that if tt is a potential function to the optimal transportation with 
positive mass distributions / and g so that u is uniquely determined up to a constant, 
the boundary point case xq G dQ can be reduced to the case xq G O by extending u to 
larger domains, and it is not necessary to define the mapping d~^u and d^u on boundary 
points. 

Taking account of our remark at the end of the previous section, we also see that 
(2.25) extends to A3w costs and uniformly c-convex domains. Consequently we obtain 
an alternate proof of the c-convexity of the solutions in Section 6 of [TW]. Futhermore 
by domain approximation, we may then extend (2.25) further to A3w costs and c-convex 
domains. 

2.6. Potential functions. Let {u,v) be potential functions to the optimal trans- 
portation problem (1.2). Then for any point xq G il, by (1.4), there exists a point yo G fi* 
such that u{xo) + v{yo) = c{xo, yo)- Hence 

(2.28) h{x) = c{x, yo) - v{yo) 
is a c-support of u at xo, and 

(2.29) h*iy) = c{xo,y)-u{xo) 

is a c-support of v at yo- 

If u is at xo, it has a unique c-support at xo- As a potential function is semi-concave, 
it is twice differential a.e. in Q. Hence u has a unique c-support almost everywhere. 

Next we consider the case when u is not at xq- Let us first introduce the terminology 
extreme point. Let i? be a convex set in M^. We say a point p G i? is an extreme point of 
E if there exists a plane P such that P HE contains only the point p. It is easy to show, 
by induction on dimensions, that any interior point in a convex set can be expressed as 
a linear combination of extreme points. 
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If u is not at an interior point xq, since u is semi-concave, d'^u{xo) is a convex set 
of dimension k for some integer 1 < A; < n. Let N:^{xo) denote the set of extreme points 
of d~^u{xo). Let 

(2.30) T^ixo) = {y e : c,{xo,y) e K{xo)}. 
Then for any yo G T^(a;o), the function 

h{x) — c{x, yo) + a 

is a global c-support of u at xq, where a is a constant such that h{xo) = u{xo). 

This assertion follows from a similar one for concave functions, which can be proved 
by blowing up the graph of w to a concave cone. That is for any p G T^(xo), there exists 
a sequence of C^-smooth points {xk} of -u, Xk — ^ xq, such that the c-support of u at 
Xk converges to a c-support of u xq. Recall that at points, u has a unique global 
c-support. 

From the above assertion and (2.7) it follows that a potential function is a local c- 
concave function if (A3) is satisfied. That is for any y' e d^u{xo), the function 

h{x) = c{x, y') - c{xq, y') -h u{xq) 

is a local c-support of u at xq. By (1.4), is a global c-support if 

(2.31) u{xQ)+v{y') = c{xo,y'). 

We remark that in general a potential function may fail to be locally c-concave (with 
respect to the definition in §2.1), if (A3w) is violated. 

From the above assertion it also follows that if u is not at xq, then the function h* 
in (2.29) is a c-support of v at any point j/o € 2^u(^o)- In other words, T^{xq) is contained 
in the contact set 

(2.32) C = {ye^*: h*{y) = v{y)}. 

By (2.31), a local c-support of u at xq is a global one if and only if the contact set C is 
c-convex (with respect to xq). 

Since Vt* is c-convex with respect to fi, from the argument in §2.5, any local c-support 
of V is a global one. Hence v is fully c-concave in Q,* . Furthermore, u is smooth if 
and only if its dual function v is strictly c-concave. 

2.7 Proof of Theorem 1. In the proof of Theorem 1 we will use the Perron lifting. 
Consider the Dirichlet problem 

(2.33) det{Dlc - D'^w) = ip(x, Dw) in Br{yo), 

w = (fi on dBriyo), 
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where ip > 0, (p^ip E C°°. Suppose r > is sufficiently small and there is a supersolution 
to (2.33) (note that when r is sufficiently small, there is always a supersolution). From 
[MTW], there is a solution w E C°^(Br{yo)) to (2.33) such that the matrix (Dlc-D'^w) > 
0. By approximation and the interior a priori estimates [MTW], there is a solution 
w e C°° {Brivo)) nC^ {Br{yo)) eC^ and ip e C^'^ Obviously w is locally c-concave. 
By §2.5, w is c-concave in B^.. 

Let (u, v) be the potential functions to (1.2). To prove Theorem 1, it suffices to prove 
that V is strictly c-concave. By approximation we assume that f,g are positive and 
smooth. By §2.4 we may also assume that Q* is uniformly c-convex with respect to Q. 
From §2.6, v is fully c-concave and every local c-support of w is a global one. Suppose to 
the contrary that v is not strictly c-concave. Then there is a c-support h* of v at some 
point yo G fl* such that the contact set 

C* = {yen*: h*{y) = v{y)} 

contains more than one point. 

From the argument in §2.5, C* cannot contain more than one disconnected component. 
In other words, C* is connected. Hence for any r > small, the intersection dB2r{yo) HC* 
is not empty. Let {Sj : j = 1, • • • , /c} be finitely many balls with radius r/2, centered on 
dBriyo), such that [j'^j^^ Bj D dBr{yo). 

Denote vq — v. For j — 1, ■ ■ ■ , k, let vj he the solution of 

(2.34) det{Dlc - D'^w) = ^So in Si, 

w = Vj-i on dBr{0), 

where 

f(x) 

(2.35) 5o = inf IdetCijI ^ ' 



which is positive by assumption of Theorem 1. Extend Vj to the whole Q,* such that 
Vj = Vj-\ in Q.* — Bj. Then v =: Vk is locally c-concave in fi*. From §2.5, it is fully 
c-concave. 

By the a priori estimates in [MTW], Vj is smooth in Bj, for j = 1, ■ ■ ■ , k. By the 
comparison principle, we have vi < vq in Bi and by induction, Vj < vq in Bj for all 
J = 2, ■ ■ • , /c. It follows that v < v near dBr{yo)- 

We have therefore obtained another fully c-concave function v which satisfies 

v = v in Br/2{yo)^{^* - Bsr/2{yo)}, 
V < V near dBr{yo). 

Hence the contact set {y E fl* : h*{y) = v{y)} cannot be connected. But this is 
impossible from the argument in §2.5 (as remarked at the end of §2.5). Hence v is 
strictly c-concave. This completes the proof of Theorem 1. 



2.8. Remarks. First we have the following result which follows from the regularity 
in Theorem 1 and approximation. 
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Corollary 1. Suppose that f,g>0, f,g& , c satisfies (Al)-(A3), and Q.* is c-convex 
with respect to Q. Then the potential function u is fully c- concave. 

Corollary 1 is a complement to the paper [MTW] . In [MTW] we introduced a notion of 
generalized solution to the boundary value problem (1.6) and proved interior regularity 
under the conditions in Corollary 1 and assuming also the smoothness of / and g. If the 
potential function u is not fully c-concave, a local c-support of u may not be a global 
one. In such case, the definition of generalized solution in [MTW] is not proper and the 
comparison principle may not hold in arbitrary sub-domains. Corollary 1 rules out the 
possibility provided the cost function satisfies (Al)-(A3) and fl* is c-convex with respect 
to f2, as assumed in [MTW]. However to get the full c-concavity of u we have to prove 
Theorem 1 first. Clearly a short and direct proof is desired. But in dimension 2, the full 
c-concavity is a direct consequence of (2.7). 

Corollary 2. Suppose the cost function c satisfies (A1)-(A3). Then any potential func- 
tions (u, v) defined in the whole are fully c-concave. 

Indeed, for any point xq & Q, if u is at xq, there is a unique global c-support of u 
at Xq. Otherwise, let yo, y\ be any two points in d'^u{xQ). Let hQ{x) = c{x, yo) — ao and 
hi{x) = c{x, yi)—ai be two c-supports of it at xq. Denote jV = {x & R" : ho{x) = hi{x)}. 
A/" is a curve which divides R^ into two parts, and both are non-compact. It suffices to 
show that 

(2.36) ht > ho (= hi) on Af, 

where t e (0, 1) and ht is as in (2.7). But if (2.36) is not true, by moving the graph of ht 
downwards, we see that there is a constant b and a point x' E Af such that 

(2.37) htix') - b = min{/io(x'), hi{x')}, 

(2.38) ht{x) — b < m.m{ho{x) , hi{x)} x & J\f near x' . 

But this is in contradiction with (2.7) at x'. Hence Corollary 2 holds. 

Note that for potential functions on bounded domains, by the uniqueness of potential 
functions when restricted to {/ > 0}, we see that u (and similarly v) is fully c-concave 
when restricted to {/ > 0}. 

The proof of Corollary 2 does not extend to higher dimensions, as we don't know 
if there is a point x' E Af such that (2.37) and (2.38) hold. But Corollary 2 holds on 
compact manifolds of any dimension, as the set J\f is compact. In particular it holds for 
the reflector design problem (in the far field case) [Wl, W2]. But for the reflector design 
problem, a c-support is a paraboloid with focus at the origin and one can also verify 
(2.36) directly [CGH]. 

We remark that if the cost function c does not satisfy (A3), then (2.7) may not hold. 
Loeper [L] shows that if condition A3w is violated, there are potential functions which 
are not fully c-concave. Furthermore, the potential function u may not be smooth 
even if both / and g are positive and smooth. Note that Loeper's potential functions are 
the negative of those here so that our c-concavity is equivalent to his c-convexity. 
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Remark. After this paper was finished, we learned that Kim and McCann found a direct 
proof of Corollary 1 above. They proved that under A3 (A3w, resp.), -^ht > (> 0, 
resp.), from which it follows that the contact set of the potential function u with its 
c-support is connected, where ht was given in (2.6). 
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